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Abstract 

To explain the anomalously large decay rate of S + — > p + it was proposed 

that a new mechanism where a light CP-odd pseudoscalar boson of m A o = 214.3 MeV 
makes a crucial contribution. Later, some authors have studied the transition 7r° — > e + e~ 
and T — > in terms of the same mechanism and their result indicates that with 
the suggested mass one cannot fit the data. This discrepancy might be caused by ex- 
perimental error of S + — > p + because there were only a few events. Whether 
the mechanism is a reasonable one motivates us to investigate the transitions 7r° — > 
e + e~; 77(77') — > i] c — > f]b — > t + t~ within the same framework. It is noted 
that for 7T° — > e + e~ , the standard model (SM) prediction is smaller than the data, whereas 
the experimental central value of 77 — > fi + fi~ is also above the SM prediction. It means 
that there should be extra contributions from other mechanisms and the contribution 
of A\ may be a possible one. Theoretically calculating the branching ratios of the con- 
cerned modes, we would check if we can obtain a universal mass for ^4^ which reconcile 
the theoretical predictions and data for all the modes. Unfortunately, we find that it is 
impossible to have such a mass with the same coupling \gi\. Therefore we conclude that 
the phenomenology does not favor such a light A®, even though a small window is still 
open. 

PACS numbers: 12.60.Jv, 13.20. Gd, 13.35.-r, 14.80.Da. 
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1 Introduction 



Searching for new physics beyond the standard model (SM) is the goal of not only 
the very high energy experiments such as at LHC and even the future ILC, but also the 
machines of lower energies where new physics signals may be revealed at rare processes. 
The HyperCP collaboration observed an anomalously large decay rate of X — > p + 
which is higher than the prediction of SM by several standard deviations [T]. The discrep- 
ancy may be attributed to contributions from new physics. A very possible mechanism is 
that the extra contribution is due to a light CP-odd pseudoscalar. Indeed, many models 
suggests its existence. Among all the models the supersymmetric model where there are 
five Higgs bosons remain after the symmetry breaking is the most favorable one. Namely 
there are two CP-even scalars H° and h°, a CP-odd A and two charged Higgs bosons 

The search for SM Higgs boson has already spanned for almost half century and 
covered a rather large regions. Recently, at LHC an excess at 126 GeV is observed and it 
could be the signal of Higgs, even though firm identification still needs time[31 HJ |5]. So 
far, as well known, the SM Higgs must be heavy, but a Higgs demanded by new physics 
beyond the SM might be light. Phenomenological search for beyond SM Higgs would be 
an interesting job of theorists and experimentalists. For example, Kao et al investigate 
the FCNC process t — » c<p with the two Higgs Doublet Model at LHC [5] where a heavy 
scalar or pseudoscalar Higgs of about 130 GeV contributes. On other aspect, it does 
not exclude the possibility that a light CP-odd pseudoscalar boson might exist, but it 
definitely is not the SM Higgs. To explain the large decay rate X — > p + which 
should be very small if only the SM applies, He, Tandean and Valencia [7] suggested that 
a light CP-odd A\ of mass 214.3 MeV may result in the observed data. Later Chang and 
Yang applies the same mechanism to evaluate the branching ratios of tt° -> e+e" while 
considering a constraint from B(T — > 7^5) [8]. They noticed that the SM prediction is 
lower than the experimental data[9], therefore there should be some extra contributions 
from the mechanisms which have not been considered yet or are due to new physics 
beyond SM. Chang and Yang calculated the contribution of the light A®, but combining 
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the constraints from the anomalous magnetic moment of muon and B(T — > 7^4?) on A®, 
the authors concluded that rigorous constraints on the mass of A\ and the concerned 
parameter \g^\ enforced by B(n° — > e + e~) and B(T — > 7^5) rule out the mass of 214.3 
MeV. Furthermore by fitting data, if the mechanism does make a substantial contribution, 
one should have m A o ~ m n and \gi\ = 0.10 ± 0.08. 

Considering that the HyperCP collaboration only recorded a few events for X — > 
p + thus relatively large experimental uncertainties could be expected, it is natural 

to ask if a mass range of A\ at vicinity of m n can remarkably enhance the rate of S — > 
p + fi + ^7 Moreover, since the new contribution to ir° — > e + e~ is realized via mediating 
a Higgs-like boson in the s-channel, its coupling to the lepton is proportional to its mass, 
so that the contribution is suppressed by the electron mass. And due to the phase space 
restriction, tt° cannot decay into heavier muons. 

In fact, we notice that for similar decay modes, the SM predictions on rj — > fi + fi~ 
is below the experimental central value, and there are no data available yet for the 
modes 7/ — » fi + fi~ r] c — > 77^ — > t + t~ etc., and we will show below that they are 

important for determining if the scenario of A\ works. The observation may hint that 
there could exist an unknown mechanism(s) which can make up the gap between the SM 
prediction and data. Existence of a light A\ definitely is a reasonable candidate. Thus 
in this work, we are going to carry out a wider study on the the modes in terms of the 
theory which involves a light CP-odd boson A\ originating from the NMSSM theory [TO] 
and the concerned couplings with fermions is: 

( - \ iA° 
£-A<iqq = ~ I / ] l u m u u>y 5 u + 22 ldm d d>y 5 d J — - , (1) 

\u-type d-type / 

£aiom = ^75^1 , (2) 

V 

where, Id = —gt = v5-/(y/2x) and l u = Id/ tan 2 ,5(7]. 

Our strategy is whether we can find a mass range as well as the parameter \gi\ (see the 
text for detail), which can tolerate all the observed modes, namely a universal A\ mass 



lr The experimental error for B(ij — > ) is large, so that making a definite conclusion needs more 

precise measurement which will be coming soon. 
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Figure 1: The QED contribution to the 77 — > £ + £ decay 



can make the gaps between the SM predictions and the data. 

The paper is arranged as follows. In section 2, we present the necessary theoretical 
derivations. In section 3, our numerical results are shown in relevant tables and figures. 
We reserve the last section for our discussion and conclusion. 

2 Formalism 

To serve our aim of this work, we concentrate ourselves on the application of the 
light CP-odd pseudoscalar boson A\ in the NMSSM. Since the leptonic decays of the 
pseudoscalar mesons tt°, rj, rj', 7] c , rjb are less contaminated by the non-perturbative 
QCD effects, they are ideal for studying the new mechanism. As aforementioned, unlike 
7T° — > e + e~, the contributions of A® to the decay modes r](rj') — > r\ c — > and 

r]b — > t + t~ do not severely suffer from the mass suppression. Then, pre-assuming the 
new mechanism, by fitting data we would check if we can obtain a universal mass for A\ 
which reconciles all the modes. 

2.1 For leptonic decays of light pseudoscalar mesons 

In the SM sector, the dominant contribution to r\ — > fi + ^ comes from the QED 
anomaly and the Feynman diagram is shown in Fig.([T]). For completeness, we re-derive 
the formulae given in Ref. [SI IUJ and show them in this text. The total contribution of 
the triangle-diagram is written as: 



M.~„ = ie 



2 



I 



d D q U w H t 



(3) 



(2tt) d (q 2 + ie)((q - p) 2 + is)((q - h) 2 - m e + ie) ' 
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where 



IT = ufasWifa-i + mtWvfas'), (4) 
= ie 2 e llUQ f 5 q a (p - qf ' f^F^q 2 \ (p - qf) . (5) 

Explicitly, is the effective 7777 vertex where the Lorentz structure includes a form 
factor related to the loop integral. The form factor, as usual, can be decomposed into a 
numerical coupling constant / 77 = l/(47r 2 / r? ) times a function F vll (q 2 , (p — q) 2 ). 

For the lepton part L^ u , we employ the projection operator technique [11] : 
V(p-ki,kx) = ■^=[v(p-ki,+)®u(k 1 ,-)+v(p-ki,-)^u(ki,+)} 

~ -^—.[-2m £ p^ 5 + \ W ( k i(P - h) T - (p - fci)V) ^ 



2vV 

+ A 5 ]- (6) 

Then we have the amplitudd^l for Fig. ([I]): 

Mjj(r] -> tt~) = -2\f2a 2 em m t m n f^A\m 2 ^) . (7) 

where A e (p 2 ) is the reduced amplitude: 

a, 2n 2 * /" A q 2 p 2 ~{q-p) 2 



A (P) = — / 2 ? 2 i \2 i ■ \?7 2 T rr -^VryCg > (P ~ ?) ) ■ ( 8 ) 

This is the same as the result given in Ref.|llj. and in the derivation, we also utilize the 
same form factor F^^q 2 , (p — q) 2 ) therein. 

Besides the QED contribution, there exists a tree level contributions induced by ex- 
changing weak interaction gauge boson Z° and a new CP-odd pseudoscalar b oson at 
s-channel and the Feynman diagrams are shown in Fig. (J2]). 

M z o = (£ + r\IV ?ez0 ~l qV, qZ oq\ V ) , (9) 
p Z * 

where q stands for the light quarks u, d and s, p = k\ + Vg ez o and Vg qZ o are the 
interaction vertices of ££Z° and qqZ° [12j . 



2 Here, we have added the missing minus according the corrections to M. E. Peskin's QFT book: 
http:/ /www. slac.stanford.edu/ mpeskin/QFT.html^crrors. 
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Figure 2: Tree level contributions to the rj — y £ + £ decay 



It is well known that the physical pseudoscalar particles rj and 7/ are mixtures of the 
flavor eigenstates r\ q and r/c[T3]: 

rjq = —r=(uu + dd), (10) 
v2 

Vs = ss , (11) 



as 



T] \ ( cos <p — sin (p \ I Vq 



(12) 



rj j \ sm<p COS</> J \ T] s 
where <ft is the mixing angle. 

Our conventions of the decay constants ffj and /* are taken from Ref . [T3| [T3] . 

(0\JU\v)=iffa 0' = ?.*). (13) 

where g stands as the lighter quarks w and d, f% = cos0/ 9 , /* = — sin0/ s and p M is the 
four-momentum of 77. 

After a straightforward calculation, we obtain the contributions of the weak interaction 
sector to the amplitude: 

2V2e 2 f^m e m r] 1 

Mz ° ~ (sin^cosMV-™! ' 
_ 2^e 2 /«m<m t? 1 

7 Z ° ~~ roinfl ^2" ' l i4 J 

(Sill Uyy COS C/p^ J p — Ul z 

where the projection operator for outgoing lepton pair is employed. 



In the NMSSM, the light CP-odd pseudoscalar Higgs A® couples to up-, down-type 
quarks and leptons. Following the general notation of Ref. [7], one can write the amplitude 
in terms of the effective couplings Eqs.flT]) and (J2]). As generally suggested in literature 



that by fitting available data, tan (3 takes a larger value, thus the coupling constant l u in 
Eq.flT]) is much suppressed and the contribution of u— type quarks to the amplitude in the 
NP part is negligible. Then the extra contribution of for d— type quarks to the amplitude 
reads: 

M A° = o AT 9 2 r— " 75^2 , (15) 

where T A o is the total width of A\. When deriving Eq. (ll5p . we utilize the relation, 

f^Tn 2 f s ui 2 
(0\d l5 d\ V ) = -i^j^ , (OMV) = , (16) 

where f% and /* are defined in Eq. ffT3j) . 
We can further reduce M A o into: 



ldgi f%m £ m 3 ri 

2v 2 m 2 - m 2 A o + im A oT A o 



M A» ~ TTTo ~2 ™2 , ~F " • l i7 J 



By the aforementioned notation, one can easily obtain A^^o as v^-M^oy?. 
The total contribution is a sum of all the individual ones: 

M tot = M, n + M u z f s + e ie ^M A f s , (18) 

where 6 NP represents a possible relative phase between the contributions of SM and 
NMSSM. 

The total decay rate of r\ — > /i + pT is expressed as: 

r to t(v -> /^ + /0 = ^-%M tot \ 2 , (19) 
Sir 

where k is the three-momentum of one of the leptons in the rest frame of 77. 

2.2 For decays of heavy pseudoscalar mesons 

Generally, when calculating the anomaly and decay rate of 7r° — > 77, for simplification, 
one can use an approximation q — > [15J. This approximation works well for decays of 
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light pseudoscalar mesons, but definitely not for heavy pseudoscalar mesons with q 2 = 
M 2 3> where M stands as the mass of the heavy meson. Thus, we take another approach 
to take into account the effects induced by the hadronic structure of the decaying heavy 
meson. 

For decay of rjb into lepton pairs, we employ the light-cone distribution amplitude 
(LCDA) method to calculate the transition amplitude M. 7J . The leading-order contribu- 
tions induced by the photon- Fermion loop are displayed in Fig. 

" — > — iWVWVWi — ' — t- ' 



WWVWv 
(-1) 




Figure 3: The QED contributions to the rjb — >■ £ + £ decay via box diagrams. 

Again, for completeness, we re-derive the QED contribution to rjt, — > t + t~, and the 
corresponding formula was obtained for i] c — > fi + fi~ in Ref.[12]. The contribution of the 
loop diagrams is: 



M 



77 



CMh^vih) + C 2 u{k l )a^v{k 2 )e, vp(7 k p 1 k' 7 



2 i 



(20) 



where 



Ci 



87T 2 J{) 



fr^me / dv4(u,fi) / dx 



l-x 



dy 



2 l + (2u-l)x-y 1 + (1 - 2u)x - y 



Co 



~^2~fvb m i J du<j){u,n) 



f l-x-y-z 

<l-l' I l dZ -jy^ 











(21) 
(22) 



Here, the notions -Di,2,3,4 and the concrete expressions for C\ and C 2 are explicitly 
presented in Ref.[12]. For 775 — > t + t~, the numerical results of C\ and C 2 are displayed 
in Table [TJ 
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Ci 


c 2 


r] b -)■ r+r 


4.22 x 10" 7 + 5.44 x 10~S 


-2.39 x 10" 8 + 2.0 x 10" 8 i 



Table 1: The coefficients C\ and C2 in Eg. ( 120]) for rjb — > t + t . 

Our next step would be evaluating the hadronic matrix element (0\b a (x)bp(y)\r]b) . 
Thus we need the wave function for the pseudo-scalars. According to Refs.[T ^ [T6 l lT7]. we 
employ the light-cone distribution amplitude for the calculation: 

(0\Kx)cMy)p\Vb(p)) = -\U b J\ue-<^ x+ ^ 0(1*, A*) , (23) 

where u = 1 — u, n is the energy scale and f Vb is the decay constant of rjb defined in 
Eq. ffTBl . The wave function of f]b is adopted as: 

cj)(u) = NAu(l - u)e~ 4 ^-^ , (24) 

where iV is the normalization factor and the parameter is set as /3 = 3.8 ± 0.7 



3 Numerical Analysis 

Firstly, we list some necessary input parameters which are taken from either the PDG 
bookp] or concerned literatures [T9l l20l [2T] |22] : 

P = 1.07/., r = 1-34A, f n = 130MeV, 

f Vt = 0.705MeV, m Vb = 9390.9 ± 2.8MeV, T Vb = lOMeV, 

= (39.9 ± 2.6(exp) ± 2.3{the))°, sin 2 6 W = 0.23116. 

The reduced amplitude A e (m^) in Eq.® is: 

A\m 2 v ) = 2.61 - 5.2U , A l {rr%) = 7.59 + 3.41i . (25) 

For the New Physics part, we firstly discuss the relevant input parameters. The elec- 
troweak scale v is 246 GeV. Furthermore, the coupling constant |^| in Eq.(j2]) is stringently 
constrained by the muon anomalous magnetic moment [7]: 

\9i\ < 1-2. (26) 
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This bound results in an A\ width T A o < 3.7 x KT 7 MeV. That implies that the coupling 
constant Id is also of order of unity which is consistent with the general estimates for the 
size of v5_/{V2x)$5\. 




1.264 x 10"' 



Figure 4: The dependance of B tot (r] — > /U + /i~) and £> iot (r/ — > on the NP phase. 

The horizontal sold lines correspond to the results of the SM, respectively. 



Generally, a relative phase between the SM and NP pieces can exist in the Lagrangian, 
and 8np should be be determined by fitting data or from a larger symmetry which includes 
the SM and the concerned NP, in this work we treat it as a free parameter. We illustrate 
the dependence of £> tot (?](?/) — > on this phase in Fig. @. From this figure one can 

observe that, as setting \gi\ = 0.35, when 6^p is around 0, Bsm(v A* + /0 is enhanced by 
the NP effect, whereas if O^p is about tt, it decreases. Comparing the result of r\ — \ fi + fi~ 
with the experimental data 

B Exp {rj^ = (5.8 ±0.8) x KT 6 . (27) 

whose central value is above the SM prediction by about 20%, we are tempted to conclude 
that as the coupling constant is not very large (\gi\ < 1.2) as suggested in the literature, 
if one expects to substantially enhance theoretical prediction of Bsuijl — > A t+ A i_ ) to the 
experimental value [18J via the effect induced by A®, the phase 6 NP should be around 
0. Thus in following calculations we set 6 NP to be 0. By contrast, Bsm(v' ~^ A t+ A t_ ) 
decreases near and almost reaches the maximum at tt. Unfortunately, up to now, there 
is no any experimental data on r/ — > , with 6 NP ps 0, its branching ratio is predicted, 
so that the future experiments may hint us if the scenario works. 
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Therefore, by fitting the experimental data of Be X p(v ~ * A t+ A i ) = (5-8 ± 0.8) x 10 6 
while taking \g^\ = 0.35 and O^p = 0, the proper mass of A\ should be 547.5 MeV. 

In order to give a better insight, we draw the dependence of B to tij]{rj') — > on 
\gt\ and m A o in Fig. (jSJ). 




1 2 5 10 20 50 100 1 2 5 10 20 50 100 200 



\gf\ \gf\ 

Figure 5: The dependance of B tot (r] — > and B to tij]' — > on the coupling 

constant |<^|. For the r] decay, the dashed line, black line and dot-dashed line correspond 
to m A o = 214.3 MeV, 500 MeV and 1 GeV respectively. For the rf decay, the black line, 
dashed line and dot-dashed line correspond to m A o = 547.5 MeV, 1 GeV, and 2 GeV 
respectively. The phase angle is chosen as 6np = 0. The abscissas are the results of the 
SM. 

From Fig.([5]), we observe that as \gi\ is not very large and 8np = 0, the NP effect for 
f] — > fi + [i~ induced by existence of A\ can enhance the branching ratio to the experimental 
level, but a heavier A\ with a mass about 1 GeV or more would not. 

Based on these parameters employed in above text, we present all the numerical results 
and corresponding experimental data for r](r]') — > [i + fi~ in Table. (J2]). 

B(r] c — > in SM was studied in Ref.fT2] and its result is: 

Bsm(Vc ->• ^ + ^) = 6.39lS;°jj x 10~ 9 . (28) 

With the same method we calculate B(f]i, — > r + r _ ) in this work. However, since Z u = 
Id/ tan 2 fl, for a larger tan/3 which is usually considered in literature, the effect induced 
by A\ on B(r] c — > /i + /x~) is negligible. Indeed, if there were the experimental data for the 
decay of r] c — > one could gain more information about A\ by comparing it with 
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Then, using m A o = 547.5 MeV, we obtain the branching ratio predicted by the pure 
SM and NMSSM with this light CP-odd Higges respectively. In analog to Fig.fll]) and 
Fig.®, we draw Fig. (jSJ) for the heavy pseudoscalar meson r\ h . Additionally, we give the 
numerical results in Table [2j 




Figure 6: (Left)The dependance of B to t(Vb — ► t + t~) on the NP phase angle. The horizon- 
tal lines represent the results of the SM. (Right)The dependance of B to t(Vb t + t~) on 
the coupling cosntant|^|, where the black line, dashed line and dot-dashed line represent 
the figures with m A o = 547.5MeV, m A o = lOGeV and m A o =30GeV, respectively. The 
phase angel is chosen as 9np = 0. The horizontal lines are the results of the SM. 

Combining with the upshots of Ref.|8J on 7r° — > e + e~, we list all the results of the 
leptonic decays of the pseudoscalar mesons 7r°, r], r/, r] c , r/b in Table. ()2]). 





&Exp 




D tot 


\m\ 


m A o (MeV) 


7T° — > e + e~ 


(7.48 ± 0.38) x 10~ 8 


(6.25 ±0.09) x 10~ 8 


7.48 x 10~ 8 


0.35 


134.95 


r] — > 


(5.8 ±0.8) x 10" 6 


4.94^ x 10- 6 


5.80 /ft x 10" 6 


0.35 


547.5 


T)' — > 




(1.27 ±0.42) x 10" 7 


(1.26 ±0.42) x 10" 7 


0.35 


547.5 






6.39^ x 10- 9 


6.39™ x 10- 9 


0.35 


547.5 




<8% 


5.56+°;^ x 10- 9 


6.671°;^ x 10- 9 


0.35 


547.5 



Table 2: The branching ratios of pseudoscalars to a leptonic pair in SM and in NMSSM. 
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jWj o(MeV) 



Figure 7: Comparison among the leptonic decays of the pseudoscalar mesons, where we 
have normalized each branching ratio with its SM prediction. The dot-dash line represents 
the experimental data and the shadowed region corresponds to the error band, where 
M = 0.35. 

4 Discussion and Conclusion 

The starting point of this work is that the SM prediction on B(tt° — > e + e~) is smaller 
than the data by a few standard deviations, and the large rate of E — > p + yU + yU~ is also 
beyond the SM prediction, thus there must be something new. But what is it, new physics 
or a mechanism hidden in the SM but not being taken into account? Definitely, it is not 
due to the final state interaction because the leptonic and semi-leptonic decays are not 
contaminated by non-perturbative QCD effects. Thus people are inclined to attribute the 
discrepancy to new physics effects. But then what new physics could it be? 

To explain the anomalously large decay rate of E + — > p + yU + /i~, He, Tandean and 
Valencia proposed a new mechanism where a light CP-odd scalar boson m A o = 214.3 
MeV exists[7]. Later, some authors studied the transition tt° — > e + e~ in terms of the 
same mechanism and their result indicates that the suggested mass cannot fit the data[8j. 
This discrepancy might be caused by experimental error of E + — > p + because there 
were only a few events. Whether the mechanism is a reasonable one motivates us to 
investigate the transitions n° — > e + e~; r\ — > r( — > r\ c — > r/ b — > t + t~ 
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within the same framework. 

Looking at Fig. (JTJ) , one can notice several aspects. 

1. Only the mass of A\ is close to the mass of the decaying pseudoscalar boson, the 
branching ratio of the leptonic decay can be remarkably enhanced. That is due to the 
Breit-Wigner form of the propagator: 

1 

q 2 -m 2 A o +imA^Al 

and q 2 = m^^^ b . 

2. Unless the mass of the light A\ is close to the mass of the decaying pseudoscalar, 
its contribution to the leptonic decay is not sensitive to the mass of A\ at all. On the 
right panel of Fig.([7j), the abscissa corresponds to the decay of r] b — > t + t~, which does 
not vary with respect to the change of m A o , even though the contribution of A\ exists. 

3. The experimental central value of B(rj — > is larger than the SM predictions, 
so that it implies that there could be additional contributions from some mechanisms 
which were not taken into account or from new physics beyond SM. On other aspect, the 
error is large, i.e. within two standard deviations, the SM prediction still coincides with 
the data. Thus if we take the central values seriously, we should search for new sources of 
the deviation. To further study, more accurate measurements are necessary. Moreover, so 
far, there are no data on r/' — > /j + available, the reason is obvious that F(j]' — > fi + ft~) 
and T(rj — > have the same order of magnitude, but T tot (r]') ^> T to t(v)- 

4. The existence of a light A\ might be the source, but our numerical results on the 
various leptonic decays of rj, rj', r/f, as well as ir° indicate that we cannot find an universal 
mass for A\ which can make the gaps between SM predictions and data for r(7r° — > e + e~) 
and T(t] — > simultaneously. 

5. Just as Chang and Yang indicated, if the mass of A\ is close to pion mass, the gap 
between theoretical prediction and data may be filled out, but another serious problem 
is raised. Namely, as 7r° and A\ have close masses, they should maximally mix according 
to the general principle of quantum mechanics, if so, the data on 7r° — > 77 would not be 
explained. Similarly, the argument can be applied to other pseudoscalar mesons. However, 
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from another aspect, as 7r° and A\ have the same CP behavior and are close in mass, in 
the two-photon final state, it is hard to distinguish between them. 

6. It is noted that the concerned experiments have larger errors, so that one may 
re-consider if we can reconcile the experimental data and theoretical predictions with the 
help of A\. Since the measurement of tt° — > e + e~ has a smaller error, let us assume that 
the mass of A\ obtained by fitting B(tt° — > e + e~) is the right one, then with this value 
we re-examine — > p + and we obtain it as 1.16 ~ 1.19 x 10~ 7 . Using the 

measured data (3.1^g ±1.5) x 1CT 8 , He, Tandean and Valencia got m A o as 214.3 MeV. 
The estimated branching ratio with a lighter A\ is 4 to 5 times larger than what they 
estimated. As pointed above, there were only a few events, the errors may be large, so that 
we hope that our experimentalists can strive to obtain more accurate measurements on 
jB(£ — y p + which may provide valuable information about A\. Moreover, if m A o is 

close to 140 MeV, as we show in Fig.([7j), its contribution to the amplitude of r\ — > ji + fi~ is 
almost a constant and negligible because the coupling of A\ to is proportional to m M 

which is small compared to m T , meanwhile the measurement on 77 — > [i + fi~ also possesses 
a larger error range and within 2 standard deviations the SM prediction is consistent with 
the data, thus that data cannot exclude an A\ of about 140 MeV. There are so far, no 
data for rj' — > fi + fi~ , 1] C — > fi + n~ and rjb — > t + t~ available, as we show, an A\ of about 140 
MeV does influence their branching ratios, but remains as a constant. For % — > t + t~, 
the result of SM+NP is about 1.2 times larger than the SM prediction. By contrast, as 
discussed in the introduction, it does not affect r\ c — > ji + \i~ at all. 

Therefore, much more accurate measurements on S — > p + and the leptonic 

decays of the pseudoscalar mesons are indeed badly needed. 

As a conclusion, the phenomenology seems not to favor the light CP-odd A®, even 
though does not exclude its existence and there exists a narrow window. 
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